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Abstract. Population studies of EGRET gamma-ray sources indicate that there is 
a distinctive population of bright sources at low galactic latitudes. The sources have 
a distribution consistent with that of young galactic objects, with a concentration 
'/^ , toward the inner spiral arms. There is a subgroup that displays strong variability 

with timescales from days to months. Following an earlier suggestion by Kaufman 
Bernado et al. (2002), we explore the possibility that these sources could be high- 
mass microquasars. Detailed models for the gamma-ray emission that include inverse 
Compton interactions of electrons in the relativistic jets and photons from all local 
\ fields (stellar UV photons, synchrotron photons, soft X-ray photons from the accre- 

. tion disk, and hard X-ray photons from a corona) are presented. We conclude that 

' microquasars are excellent candidates for the parent population of the subgroup of 

, variable low-latitude EGRET sources. 
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1. Introduction 



Population studies of EGRET gamma-ray sources indicate that there 
is a distinctive population of bright sources at low galactic latitudes 
5^ ■ (Gehrels et al. 2000; Grenier 2001, 2004; Romero 2001). They are 

well-correlated with star forming regions and HII regions, which is 
indicative of an association with young stellar objects (Romero et al. 
1999, Romero 2001). Log - logS" studies suggest that they are more 
abundant toward the inner spiral arms (Gehrels et al. 2000, Bhat- 
tacharya et al. 2003). These are bright sources (isotropic luminosities 
in the range 10^^"'^^ erg s~^), with an average photon spectral index 
(r) = 2.18 (Gehrels et al. 2000), where F{E) oc E''^ . These sources, 
whose number is ~ 45 it 9, form the Gamma-Ray Population I (GRP 
I, Romero et al. 2004). Among GRP I sources there is a subgroup that 
displays significant variability on timescales of weeks to months (Torres 
et al. 2001, Nolan et al. 2003). Recently, Kaufman Bernado et al. (2002) 
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and Romero et al. (2004) have suggested that this subgroup of GRP 
I sources might be high-mass microquasars (i.e. microquasars formed 
by a compact object and an early- type stellar companion), where the 
gamma-ray emission arises from interactions between relativistic parti- 
cles in the jet and external photon fields, most notably the stellar UV 
emission. 

In the present paper we explore in more detail this hypothesis, 
presenting more realistic models for the gamma-ray emission. In par- 
ticular, we will include effects of the interaction of the microquasar 
jet with the X-ray fields produced by the accretion disk and the hot 
corona that is thought to surround the compact object. We will also 
include synchrotron self-Compton emission, Klein-Nishina effects, and 
the back-reaction of the different losses in the particle spectrum of 
the jet. We will calculate, for some representative sets of the param- 
eters that characterize high-mass microquasars, the spectral energy 
distribution (SED) from radio wavelengths up to GeV gamma-rays. 



2. GRP I sources 

GRP I sources concentrate along the galactic plane and present a good 
spatial correlation with young stellar objects (Romero et al. 1999). 
The variability analysis of these sources by Torres et al. (2001) clearly 
shows evidence for the existence of a subgroup with variable emission 
on timescales from weeks to months. This is corroborated by the recent 
results presented by Nolan et al. (2003), which are based on a maximum 
likelihood re-analysis of the EGRET data. These authors identify 17 
variable sources within 6° from the galactic plane. These sources are 
clumped within 55° of the galactic center. 

A log N — log S analysis for all GRP I sources yields a distribu- 
tion that is consistent with a power- law with index {(3) = 3.1 it 0.4 
(Bhattacharya et al. 2003). This is far steeper than what is expected 
for a population uniformly distributed along the galactic disk. For in- 
stance, for pulsars detected at 400 MHz the slope results ~ 1.05 
(Bhattacharya et al. 2003) . The unidentified gamma-ray sources, on the 
contrary, seem to be concentrated mainly in the inner spiral arms. In 
order to find possible further evidence for different populations among 
GRP I sources, we have implemented a log N — log S analysis of both 
variable and non- variable low-latitude sources. 

First we have considered the 17 variable sources listed by Nolan et al. 
(2003). In order to take into account systematic effects introduced by 
different exposure and background resulting in non-uniform detectabil- 
ity, we have adopted the procedure described by Reimer (2001). The 
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Figure 2. Log N — log 5 plot for non- variable gamma-ray sources within |6| < 6°. 

obtained log N — log S plot is shown in Fig. 1. The normalized distribu- 
tion can be fitted by a power-law N{S) oc S~/^, with {/3) = 1.66 ± 0.31, 
significantly flatter than for the entire sample. If we now consider those 
sources that classified as non- variable or dubious cases, we get the 
log N — log S plot shown in Fig. 2. In this case the distribution can 
be fitted by a power-law with index = 2.92 it 0.36. Quoted errors 
for P are 1 sigma errors. 
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The average spectral index is also different for both samples: in the 
case of the variable sources we have (F) = 2.04 ib 0.1, whereas for the 
remaining sources is (F) = 2.16 it 0.2. However, the averaged errors are 
still large, due to the low number of sources, and do not allow yet to 
distinguish both populations properly. In any case, log TV — log S analy- 
sis suggest that there arc two different groups of sources, one formed by 
steady sources concentrated toward the inner spiral arms, and a second 
group with variable sources and a wider distribution along the galactic 
plane, although not as wide as that of radio pulsars (Bhattacharya et al. 
2003). 

High-mass microquasars appear to be good candidates for compact 
and variable sources in the galactic plane. Since they can have very 
large proper motions (e.g. Ribo ct al. 2002), their distribution along the 
plane should be broader than that presented by supernova remnants 
and molecular clouds, which can be traced by star- forming regions and 
OB associations. Their spread, however, is limited by the lifespan of 
the companion massive star, and hence it is not as extended as that 
of radio pulsars. In the next section we will discuss the potential of 
microquasars as gamma-ray sources. 



3. Gamma-ray emission from microquasars 

In the present paper we will consider a relativistic inhomogeneous jet 
flow moving along the z axis (assumed to be perpendicular to the disk) 
with a bulk Lorentz factor Fjet- We shall allow the jet to expand later- 
ally, in such a way that the radius i? at a distance z from the compact 
object will be given by R{z) = (^z^ , with e < 1 and zq < z < 2;max- C 
is related to the opening angle of the jet. For e = 1 we have a conical 
jet. The electron energy distribution {N{z,'^e), a power-law of index 
p), the maximum electron Lorentz factor (7emax(-z)) and the magnetic 
field {B{z)) will be parametrized like in Ghisellini et al. (1985). 

The synchrotron radiation density within the jet can be estimated 
using the local approximation to the synchrotron radiation field (e.g. 
Ghisellini et al. 1985). We have used the formulae presented in Pachol- 
czyk (1970), adapted to an expanding jet. The total radiation field to 
which the leptons are exposed in the jet will have also contribution from 
external sources. These contributions can be modeled as two black body 
components, one peaked at UV energies (the companion star field) and 
other at energies kT ~ 1 keV (the inner accretion disk field), plus a 
power-law with an exponential cutoff at kT ~ 150 keV (the corona). 
With the exception of the disk, these contributions are assumed to be 
isotropic (Romero et al. 2002, Georganopoulos et al. 2002). 
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Once N{z, 7e) and the total energy density in seed photons (C/tot(€0) z)) 
have been determined, the inverse Compton (IC) interaction between 
tliem can be studied calculating the IC spectral energy distribution per 
energy unit. The cross section ((7(x, cq, 7e)) for both the Thomson and 
the Klein-Nishina regimes has been taken from Blumenthal &: Gould 
(1970). The spectral energy distribution for the optically thin case in 
the jet's reference frame is: 

(z) /-TemaxCz) 

e^e = e / / / S(z)C/tot(eo,^;) 

■2-min ^Omin Tcmin (^) 

xN{'ye,z) djdeodz, (1) 

ae eg 

where €o and e are the energies of the incoming and the outgoing 
photons, respectively. The parameter x is actually a function which 
depends on both of the former quantities and on the electron Lorentz 
factor. T,{z) is the surface of a perpendicular jet slice located at z. 
Notice that the external fields contributing to Utot{^o,z) should be 
transformed to the co-moving frame. Detailed expressions for such 
transformations are given by Dermer &; Schlickeiser (2002). 
In the observer's reference frame we have: 

e'L'^, = D'^+Pe'L,,. (2) 

The integration is performed in the co-moving system and then the 
result is transformed to the observer's frame, hence the factor D^'^^, 
which is the Doppler boosting for a continuous jet. The energy of the 
scattered photons in the jet's reference frame (e) is boosted to e' = De. 
The Doppler factor D for the approaching jet is given by 

^ " rjet(l-/?COS^)' 

where /3 is the velocity of the jet in speed of light units and 9 is the 
angle between the jet and the line of sight. 

In the case of the IC interactions with disk photons, a factor (1 — 
cos 0)(P+^)/2 must be introduced in Eq.( 1) in order to take into account 
the fact that the photons come from 6e/imdthe jet (Dermer et al. 1992). 

In order to make any calculation of the IC emission of a given mi- 
croquasar, we have to specify first the jet power. In this work we shall 
adopt the disk/jet coupling hypothesis formulated by Falcke & Bier- 
mann (1995, 1999), i.e. the jet power is proportional to the accretion 
rate: 

-^jet = QL^c- (4) 
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Here, Lac = is the accretion power onto the compact object and q 
is a number < 1. For a leptonic jet, as in our case, q <^1. The jet power 
determines the energy content of the electron energy distribution for a 
leptonic jet (see, e.g., Georganopoulos et al. 2002). 

Although the model outlined above concerns primarily to the high 
energy emission range since it focuses on the IC radiation from the 
inner jet, the synchrotron formulae can be used to make predictions 
also at radio wavelengths. With this aim, we have developed an outer 
jet model appropriate for regions farther from the compact object. This 
outer jet model will fulfill most of the conditions of the former one, but 
it will require a different geometry because of the expected lower energy 
loss rate of the particles. 

For further details about our model, see Bosch-Ramon et al. (2004). 



4. Specific models 

In order to explore a plethora of situations, we have adopted a set of 
assumptions, which are summarized in Tables I. In Table II we list those 
parameters that were allowed to change along with the values adopted 
for three representative 'mild' microquasar, a microblazar, and 

a system similar to LS 5039 and LS I +61 303^. 

In Fig. 3, we show the whole SED of the 'mild' case. We plot 
here the spectra for all contributions to the emission, i.e. the outer 
jet, the seed photon sources, and the IC components (for outer jet 
parameters, see Table III). In Fig. 4 we show a more extreme case, 
for a source with a high bulk Lorentz factor rjet=10, a high-energy 
cutoff for electrons of 7emax(-2o) = 10^, and a small viewing angle 
9 = 1° (see Table II). Finally, we plot a SED similar (realistic case 
in Tables II and III) to the one presented by two microquasars that 
arc considered as the best candidates to counterparts of unidentified 
EGRET sources, LS 5039 and LSI +61°303 (3EG 1824+1514, Paredes 
et al. 2000, and 3EG 0241+6103, Kniffen et al. 1997, respectively - see 
Fig. 5); the selected set of parameter is consistent with those adopted 
in Bosch-Ramon & Paredes (2004a, b). 



^ LS 5039 and LSI +61° 303 present significant diflPerences though, regarding 
gamma-rays, they are pretty similar. Also, at lower energies, both systems can be 
treated in a similar way at this stage. 
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Parameter (symbol) Value 



Black hole mass (Mbh) lOM© 

Gravitational radius (-Rg) 1.48 x 10® cm 

Accretion luminosity (Lac) 10~* Mqc^ year~^ 

Stellar radius {R*) 15 Rq 

Stellar bolometric luminosity (L«) 5x10^* erg s~^ 

Distance from jet's apex to the compact object (20) 50 -Rg 

Initial jet radius (Tio) O.l^o 

Orbital radius (Rorh) 45 Rq 

Peak energy of the disk (fcTdisk) 1 keV 

Peak energy of the corona (fcTcor) 150 keV 

Peak energy of the star (fcTstar) 10 eV 

Expansion coefficient of the inner jet (e) 1 

Electron power-law index (p) 2 

Photon index for the corona (Fcor) 1-6 



5. Discussion 

In a previous work (Bosch-Ramon et al. 2004), it was shown that for 
several models with disk/jet coupling constant q = 10"^ we can get the 

Table II. Changing parameters in the different cases considered in the text. 



Parameter (symbol) Value 

mild/oxtrcnic/rcalistic 



Corona luminosity (Lcor) 

Total disk luminosity (Ldisk) 

Inner jet magnetic field {B{zo)) 

Maximum electron Lorentz factor (7emax(2o)) 

Bulk Lorentz factor of the jet (Fjet) 

Viewing angle to jet's axis 6 

Disk/jet coupling constant (q) 



10^^ 10^^, 3 X 10^^ erg s"^ 
10^^, 10^'^, 5x10=^2 erg s"^ 

10, 10, 200 G 

10", 10^ 10" 

2.5, 10 , 1.1 

10°, 1°, 10° 
10-*, 10"*, 10"* 
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Figure 3. Complete SED for a model with q = 10"", B=W G, Ldisk = lO^'^ erg s~\ 
-Lcor = 10^^ erg s~^, jeinax{zo) = 10*, rjet=2.5, and a viewing angle of 10°. The 
synchrotron seed photons come from the inner jet. 




Figure 4. SED for a model with q = 10"", S=10 G, Ldisk = lO^'' erg 
icor ~ 10^^ erg s"^, 7cmax(2:()) — 10®, rjct=10, and a viewing angle of 1°. The 
synchrotron seed photons come from the inner jet. 
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Parameter (symbol) 



Value 

mild /extreme /realistic 



Expansion coefficient of the outer jet (e) 

Initial point of the outer jet (20) 

Initial radius of the outer jet (i?o) 

(Computed) maximum length of the jet (few GHz) 

Maximum electron lorentz factor of the outer jet (7ei 

Outer jet magnetic field {B{zo)) 

Electron power-law index (p) 



0.2, 0.3, 0.2 
500 i?g, 5000 i?g , 500 Rg 
50 500 7?g, 50 



ll6 



lo^ 10'', 10^ 

1, 0.1, 20 G 
2 




Figure 5. SED for a model with q = 10"^, S=200 G, Lcor and 
i-disk = 3 X 10^^ erg s~^, 7emax(^;o) = 10", rjet=l.l, and a viewing angle of 10°. 
The synchrotron seed photons come from the inner jet. 



expected luminosities in the observer's frame inferred for GRP I sources 
with the right photon index at energies ~ 1 GeV, i.e. F ^ 2. When the 
magnetic field is strong enough {B(zq) ^ 100 G), SSC emission alone 
can account for luminosities ~ 10'^^ erg s~^ at 1 GeV. Models with bulk 
Lorentz factors Fjet = 1.5 and 5 do not produce dramatically different 
results for a jet with a viewing angle of ~ 10° except for the case 
of the scattered stellar photons. In case of lower magnetic fields, the 
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IC scattering upon external fields clearly dominates. It is interesting 
to notice that even models with light jets {q = 10^^) can produce 
significant gamma-ray sources {rsj IQ-^"^ erg s^^ at 100 MeV) when a 
strong corona is present. In these models, the spectrum tends to be a 
bit harder than in the case of stellar photons, with our current set of 
assumptions (value of Tcor, etc). 

In the present study of concrete cases, we compare Figs. 3, 4 and 5. 
Prom these figures, it is possible to see the diflferences between 'mild' 
microquasars, with Tjet = 2.5, a viewing angle oi9 = 10° and q = 10~^, 
and more 'extreme' cases with Fjct = 10 and 9 = 1°. Since we have 
assumed in both cases a magnetic field of B{zo) = 10 G, external IC 
scattering always dominates over SSC. In the 'mild' microquasar the 
emission at MeV-GeV energies is mainly due to upscattering of coronal 
photons whereas, in the extreme case, the IC emission in the stellar 
field exceeds by far the other contributions. For the case similar to 
LS 5039/LSI +61°303, we have taken the magnetic field to be high 
enough to be the main contribution at high energies to reproduce the 
observed features of the spectra; mainly: higher luminosities at gamma- 
rays than at X-rays and no strong disk/corona hints at X-rays. In the 
three cases the radio counterpart is rather weaJi, with a Lradio/-f'7 ^ 1- 
The hard X-ray counterpart, on the contrary, can be important, even 
beyond the cutoff for the corona (at ~ 150 keV). Hence, INTEGRAL 
observations can be useful to unmask candidates that are obscured at 
optical wavelength (see Combi et al. 2004 for a recent study in this 
direction) . 

An interesting feature of the extreme case is that the gamma-ray 
emission with a hard spectrum F ~ 2 extends up to --^100 GeV. Above 
10^^ eV the spectrum becomes softer due to the Klein-Nishina ef- 
fect. This type of sources should be detectable with modern Imaging 
Cherenkov Telescopes like HESS and MAGIC. Even systems with low- 
mass companions and weak coronas might be detectable since the SSC 
component has luminosities of ~ 10^^ erg s~^. The soft spectrum can 
be used to differentiate this leptonic model from other hadronic alter- 
natives which predict a harder spectrum at very high energies (e.g. 
Romero et al. 2003). 



6. Conclusions 

We have shown that the variable gamma-ray sources found on the 
galactic plane have some common features that make reasonable to 
consider them as a distinctive group of GRP I sources. We have sug- 
gested that these sources might be microquasars with high-mass stellar 
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companions and we have elaborated some detailed models to explain 
the gamma-ray production in this type of objects. In particular, we 
have considered inhomogeneous jet models where gamma-rays are the 
result of inverse Compton interactions of leptons in the inner jet with 
locally produced synchrotron photons as well as external photon fields. 
We have calculated the emission resulting from the upscattering of 
disk, coronal, and stellar photons, incorporating a full Klein-Nishina 
calculation and the effect of losses in the particle spectrum. We have 
shown that different spectral energy distributions at high energies can 
be obtained from different and reasonable combinations of the physical 
parameters like magnetic field, jet power, coronal and disk luminosities, 
etc. The model can compute with a reasonable set of parameters a 
realistic SED, reproducing the ones associated to known microquasars. 
It seems to be clear that the microquasar phenomenon can be naturally 
extended up to the highest energies and that we can expect these 
objects to manifest themselves as a distinctive group of gamma-ray 
sources that might be detectable with satellite-borne instruments like 
those to be carried by AGILE and GLAST and even by ground-based 
Cherenkov telescopes hke HESS and MAGIC. 

Microquasars are not, of course, the only kind of galactic object that 
might display variable gamma-ray emission. Other alternatives include 
early-type binaries (Benaglia & Romero 2003), accreting neutron stars 
(Romero et al. 2001), pulsar wind nebulae (Roberts et al. 2002) and 
exotic objects (Punsly et al. 2000). However, microquasars arc perhaps 
the most attractive candidates to explain a significant fraction of the 
variable GRP I sources because of the well-established presence of 
relativistic jets in these objects, along with several external photon 
fields. 
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